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ABSTRACT 

We study the bar-driven dynamics in the inner part of the Milky Way by using invariant mani- 
folds. This theory has been successfully applied to describe the morphology and kinematics of 
O ■ rings and spirals in external galaxies, and now, for the first time, we apply it to the Milky Way. 

In particular, we compute the orbits confined by the invariant manifolds of the unstable peri- 
c/^ . odic orbits located at the ends of the bar. We start by discussing whether the COBE/DIRBE 

, ■ bar and the Long bar compose a single bar or two independent bars and perform a number of 

comparisons which, taken together, argue strongly in favour of the former. More specifically, 
we favour the possibility that the so-called COBE/DIRBE bar is the boxy/peanut bulge of a 
^ ■ bar whose outer thin parts are the so-called Long bar. This possibility is in good agreement 

I both with observations of external galaxies, with orbital structure theory and with simulations. 

■ We then analyse in detail the morphology and kinematics given by five representative Galac- 

tic potentials. Two of these have a Ferrers bar, two have a quadrupole bar and the last one 
a composite bar. We first consider only the COBE/DIRBE bar and then extend it to include 
00 ■ the effect of the Long bar. We find that The large-scale structure given by the manifolds de- 

I scribes an inner ring, whose size is similar to the near and far 3-kpc arm, and an outer ring, 

whose properties resemble those of the Galactic Molecular Ring. We also analyse the kine- 
matics of these two structures, under the different galactic potentials, and find they reproduce 
the relevant over-densities found in the galactic longitude - velocity CO diagram. Finally, we 
. - consider for what model parameters, the global morphology of the manifolds may reproduce 

/\ ' the two outer spiral arms. We conclude that this would necessitate either more massive and 

more rapidly rotating bars, or including in the potential an extra component describing the 
spiral arms. 
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1 INTRODUCTION 

The large-scale structure of the Milky Way (hereafter, MW) 
disc has been under study for many years. T he COBE/DIRBE 
('Weiland et al \ \l994) and Spitzer/GLIMPSE (IChurchwell et al. ] 
>2QQ9) missions provided infrared information on the global struc- 
ture of the inner Galaxy. Even though these studies have pro- 
vided some light, the large-scale structure of the MW disc proves 
to be highly complex. Near- and mid-IR low resoluti on images 
detected the COBE/DIRBE bar dWeiland et al] Il994l) . also re- 
ferred to as the triaxial bulge or the COBE/DIRBE bar. Near- 
IR red clump giants of the mid plane r evealed the existence of 
a second bar dHammerslev a/. I [200Q| : iLopez-Corredoira et al. I 
I2OO7I: ICabrera-Lavers et al. WlOOA and confirmed by GLIMPSE 
teeniamin et al. Il2005h , usuallv referred to as the Long bar. Ob- 



servations suggest two other large scale structures towards the in- 
ner parts of the Milky Way, namely the near and far 3-kpc arms 
teenlll964lDame& Thaddeu]l2008h and the Galactic Molecular 
Ring Jciemens, Sanders & Scovill3 19881) . Although their charac- 
teristics or even their existence are currently being under debate 
toame & Tha ddeus 20 11), here we aim to bring some light given 
the observed characteristics up-to-date. The near and far 3-kpc 
arms were detected using the HI 21 -cm line and CO emission sur- 
veys and extend roughly parallel to the COBE/DIRBE bar, whereas 
the position of the Gal actic Molecular Ring (hereafte r, GMR) is 
not so well determined. Iciemens, Sanders & Scovilld 11988.) sug- 
gested it is located at ^ 5.5 kpc from the Galactic Center, while 
other authors suggest it is located about halfway to the Galactic 
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Center Jiinnev & Memfieldl Il998l : IPame. Hartmann & ThaddeusI 
l200ll : lRathborne etal Il2009h . 

There has been a lot of effort to determine the bar and 
spiral arm characteristic s, both from the observational a n d the- 
oreti cal point of vie w. IPebattista. Gerhard & Sevenstej (l2002h 
and ISevenstej (l2002h use a sample of - 250 bright OH/IR 
stars of the inner Galaxy to determine the pattern speed of 
the COBE/DIRBE bar, while the CO emission and the 21 -cm 
line of n eutral hydrogen map the galactic longitude-velocity (l,v) 
diagr am ( Dame. Hartmann & ThaddeusI I2OOII : IPame & ThaddeusI 
I2OO8I : IValleell2008h . The spiral arms and rings of the Galaxy ap- 
pear as over-densities in such diagrams. From the theoretical point 
of view, several works use hydrodynamics simulations to constrain 
the Galaxy parameters and to reproduce the observed (l,v) diagram 



the ( jalaxy parameters and to reproduce the observed (l,v) diagram 
(e.g. Englmaier & Gerhardll 19991: Rodriguez-Fernandez & CombesI 



l2008:'Baba. Saito h & Wadall201Ql) . From the orbits point of view. 
^Habing^ra/. (20# compute a library of orbits, some of which 
reproduce the over-densities in th e inner longitudes of the (/,^;)- 
diagram, while iGreen et al. I (l201lh study in detail the shape of the 
3-kpc arm concluding that it can be approximated by an elliptical 
ring. Test particle simulations have recently been used to study the 
velocity distribution function in the Solar neighbourhood and to 
use it to constrain the characteristics of the bar and/or spiral arms 
(lDehnenll2000a[: BI2OOII : IChakrabartvll2007l : lAntoia g/a/.1l2009l : 
iMinchev et al iboiol) . 

In this paper we will use an approach which is novel for 
our Galaxy, namely that of invariant manifolds. In a previous 
set of papers we developed the method and techniques that 
we will use here and also applied them to th e study of spi- 
rals a nd rings in external galaxies (see iRoine ro-Gom ez a/ I 
20061. I2OO7I: lAthanassoula. Romero-Gomez & Masdemon^ l2009l : 



Athanassoula et al. I2OO9I . I2OIOI . lereafter Papers I-V, respectively). 
The invariant manifolds are linked to the presence of the La- 
grangian points Li and L2 of a barred system and they can re- 
produce the observed structures of rings and spiral arms. In pa- 
pers I-V, we studied in detail the characteristics of the orbits con- 
fined by such invariant manifolds and we analysed both their mor- 
phology and kinematics so as to compare them to the rings and 
spiral arms in external barred galaxies. Invariant manifolds, al- 
beit in a quite different way th an what we have her e and Pa- 
pers I-V, have also been used by iTsoutsis et al. 1 (I2OO9I) to model 
thr ee barred galaxies, NGC 3992, NGC 107 3, and NGC 1398 and 
by IPatsis. Kalapotharakos & Grosboll (I2OIC1I) to perform an orbital 
analysis of NGC 1300. 

Here we use the invariant manifolds to compute a family of 
orbits of a wide range of energies and to study their morphology 
and their kinematics. We want to evaluate whether the invariant 
manifolds can provide an alternative, plausible model for the inner 
part of the MW disc. By using MW analytical potentials recently 
used in the literature, we try to answer the following questions: 
Can observations be plausibly interpreted by manifolds? Can the 
latter provide an alternative interpretation of the inner structure of 
the MW? Which are the requirements the potential has to fulfill in 
order to reproduce the rings and spiral arms of the MW? In a future 
paper we will examine whether the combination of the manifolds 
and observations can constrain the MW bar properties. 

In order to compute the orbits confined by the manifolds we 
must first fix the Galactic potential. We have chosen the most rep- 
resentative analytical potential s used to desc ribe a COBE/DIRBE 
bar, namely a Ferrers b ar (iFerrersI Il877h . a quadrupole bar 
( Binnev & Tremain^ l2008h and a composite bar JPichardo et al. I 
I2OO4I) . Deliberately, we consider studies where the authors have 



tuned one of these three types of potentials to the COBE/DIRBE 
bar, and we use the same set of parameters as they do. 

The paper is organised as follows: in Sect. [2l we discuss 
whether our Galaxy has a single or a double bar, using arguments 
from the morphology of external galaxies, from orbital structure 
theory and from A/^-body simulations. We also describe the models 
and compare them in detail in terms of forces. In Sect. (3] we give a 
brief summary of the dynamics driven by the unstable Lagrangian 
points and, in particular, the definition of the invariant manifolds. 
We also give a brief summary of the main relevant results found 
in papers I-V. In Sect. IH we compute the invariant manifolds for 
the selected models and we analyse them in terms of morphol- 
ogy and kinematics. The results are compared to the observables 
in Sect. 14.31 In Sect. [S] we explore the parameter space and de- 
termine in which cases the manifolds could reproduce outer spiral 
arms. Finally, we give a short summary and conclusions in Sect. (6] 
In the Appendix, we describe in detail the analytical models and 
give the default parameters used. 



2 MODELLING THE GALACTIC POTENTIAL 
2.1 Analytical models 

There are several analytical models in the literature used to model 
the MW Galaxy. They essentially consist of an axisymmetric plus 
a one-bar component. Each model has been constructed to model 
the Milky Way and in the Appendix of this paper we give a brief 
description of the potentials used and their default parameters. We 
want to stress here that we will consider the same parameters as 
these studies. The axisymmetric component describes the disc, halo 
and bulge of the Galaxy and in each model it is modelled in a dif- 
ferent way. 

The models considered in this paper are: 

• Melnik & Rautiainen (2009, hereafter MR09) and Gardner 
& Fly nn (2010, hereafter GFIO). Both use a Ferrers bar (iFerrersI 
Il877h . though the purpose of each of the papers is very differ- 
ent. The former uses the bar potential in test particle simulations 
to model the kinematics of the outer rings and spirals of the Galaxy 
and to compare it with the residual velocities of OB -associations 
in the Perseus and Sagittarius regions. The latter studies the effect 
of the bar parameters on the kinematic substructures found in the 
velo city plan e of the Solar neighbourho od. 

• iDehnenI feOOOal) and |f3 (1200 ll) (hereafter, DehnenOO and 
FuxOl, respectively). Both use a quadrupole bar, but with differ- 
ent model parameters, to study the effect of the COBE/DIRBE bar 
on the local disc stellar kinematics. Such potentials are also often 
referred to as ad-hoc potentials, since they are not obtained from 
any density distribution. 

• The composite bar of Pichardo et al. (2004, hereafter PMM04) 
consists of a set of prolate ellipsoids, superposed so that the 
surface density matches the mass distribution obtained by the 
COBE/DIRBE mission, and from which the potential and forces 
are derived. The authors compute in detail the families of orbits 
given by this potential and use surfaces of section to characterise 
the bar structure. 

In order to fix the parameters of each model, the authors 
take into account the available relevant observational data. Even 
though there is a lot of uncertainty in these data, a range of 
possible values can be determined. The semi-major axis of the 
COBE/DIRBE bar or Galactic bulge is - 3.1 - 3.5 /cpc and its 
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aspect ratio of 10 : 4:3 (lengt h : width : height), as es ti mated 
by COBE/DI RBE dWeiland a/H 1 19941 : iFreudenreichI Il998l : 
lGerhardl[2002h . Several studies fix the m ass of the COBE/DIRBE 
bar in the range 1 - 2 x 1 ^° (iMatsumoto a/. I Il982l: 
Ki^ 1 19921 : iDwek et al 1 1 19951 : IZhad Il996l : IWeiner & SellwoodI 



19991) . The relative orientation of the bar with respect to the 
Galactic Center - Sun is not well established, although most 
observations (2MASS star counts or red clump giants) and models 
(based on COBE/DIRBE and Spitzer/GLI MPSE) agree it lie s 
roughly within the range o f 15° — 3 0° (e.^. iDwek g^ (1 19951); 
iBinn ev, Gerhard & Spergel^ 
(|1999): Endmaier & Gerhard 



f anekg^ aTI <1997b; IfuxI 
issantz & Gerhard (l2002l): 



iBissantz. Englmaier & Gerh^ ( 2003 ): Babusiaux & Gilmore' 

(l2005h: I Lopez-Corredoira. Cabrera-Lavers & Gerhard ( 2005^) ; 
iBeniamin g^ g/. I ( l2005l) : IChurchwell g^ a/. I ( 120091) ). Here we 
choose 20° as a representative value. The pattern speed obtained 
by the studies mentioned in the Introduction lies within the 
range Q^h — 35 — 60A;m g ~^ kpc~^ although higher valuers are 
favored (see lGerhardll20ld for a review). As for the Long bar, it is 
somewhat longer than the COBE/DIRBE bar with a semi-major 
axis of ^ 4 — 4.5 kpc an d an aspect ratio of 10 : 1.54 : 0.26 
jHammerslev et al I I2OOOI) . The mass of the Long bar is less 
than the mass of the COBE/DIRBE bar. It is estimated to be 
around 6 x 10^ Mq^ i.e. about 2/3 that of the COB E/DIRBE bar 
JHammerslev et a/. I l200d : [Gardner & Flvmil I2OIOI) . The relative 
orientation from the Galact ic Center - Sun l ine is estimate d obser- 
vationally to be ^ 40° (»Hamniersley g^ al. [2OOO: Benjamin et al. I 



l2005l : lL6pez-Corredoira et al. II2007I : ICabrera-Lavers et al. II2OO8I) . 
Note that the observational sources used to constrain the pa- 
rameters of the COBE/DIRBE and the Long bar are different. 
These could, in principle, lead to somewhat different definitions 
of bar length. However, according to tFreudenreich (1998) for the 
COBE/DIRBE bar and iHammerslev g^ a/. I ilOoS) for the Long 
bar, both give estimates of the bar half length. Even though these 
estimates are subject to observational errors, the ratio of the two 
values will not change significantly. 

2.2 How many bars does our Galaxy have? 

A considerable fraction of external barred galaxies are known to 
have two bars: a primary, or main bar and a secondary, or inner 
bar. This fraction depends on the galaxy's Hubble type, the qual- 
ity of the sample images, and other factors, but fractions of the 
order of a third or a fourth are quite reasonable. Could it be that 
our Galaxy is one of these? In order to pursue this line further, one 
needs to make sure that the properties of the COBE/DIRBE bar 
and the Long bar are compatible with those of galaxies with double 
bars. Can we safely assume that the Long bar is the main bar and 
that the COBE/DIRBE bar is the secondary bar? 

Several propertie s of inner bars have been well s tudied (e.g. 
lErwin & SpaSbool Eaine et al.|[2002l . lErwin|[200^ and some 
major trends have been found. The strongest constraints come from 
the bar length. Inner bars are quite small, with a semi-major axis 
between 100 pc and 1.2 kpc, with median size around 500 pc. 
Typically their rela tive length is about 12% of that of the main 
bar. The sample of lErwinI (l2009l) contains 64 galaxies with dou- 
ble bars and out of these only two have a secondary bar longer 
than 22% of the primary, and none longer than 30%. This con- 
trasts strongly with the numbers for our Galaxy, where the length 
of the bar semi-major axis are 3.1 - 3.5 kpc and 4 kpc for the 
COBE/DIRBE and for the Long bar, respectively. Thus the length 
of the COBE/DIRBE bar is more than 10 sigmas beyond the val- 




Figure 1. Two edge-on views of a bar from a simulation (see text). The 
upper panel shows a side-on view, and the bottom one a view from an angle 
near the bar major axis. In both cases the views are perspective and thus the 
relative thickness of the inner and outer parts should be inferred only from 
the upper panel. 




Figure 2. Face-on view of a simulation. This is not a model built specifically 
to represent our Galaxy, but is a clear example of a snapshot with a short 
leading ring segment emanating from the end of the bar (see text). In this 
figure the rotation is clockwise to facilitate comparison with our Galaxy. 



ues found for the inner bars of external galaxies. Furthermore, the 
relative length of the COBE/DIRBE bar relative to the Long bar is 
0.8, again more than 10 sigmas out of the distribution found from 
external galaxies. Double bars have also been found in simulations 
(iHeller. Shlosman & Athanassoulall2007l : IShen & Debattistall2009l) 
and their parameters are in good agreements with those of observed 
double bars and in disagreement with the values for our Galaxy. 

So one can reasonably exclude that the Long and the 
COBE/DIRBE bar form a double bar system, because their prop- 
erties are very far from those of double bars in external galaxies 
and simulations, making the two incompatibl^B- So then what is it? 
The data give us an important clue for that. Namely the ratio of the 
major- to z- semi-axis of the bar is 0.3 for the COBE/DIRBE 
bar and ^ 0.026 for the Long bar, i.e. the Long bar is very thin 
and the COBE/DIRBE bar is very thick. Let us therefore exam- 
ine the alternative that there is only one bar in the galaxy and 
that the COBE/DIRBE bar is simply the part which corresponds 
to the boxy/peanut bulge and the Long bar the outer part of this 
bar. This g eometry has been already discussed for external g alax- 
ies in (e.g. Athanassoula 2005; Athanassoula & Beaton 2006; and 
was first proposed for our G alaxy by Athanasso ula (2006, 2008b . 
ICabrera-Lavers et~al. ](|2003) tested this suggestion using their red- 
clump giants measurements. We will discuss here further what the 
rele vant orbital structure stu dies and iV-body simul ations imply. 
IPfenniged (Il984h and ISkokos et al \ (l2002al lbl) studied the 



1 lAlardI ( l200lh , however, found evidence for a small lopsided bar in our 
Galaxy, whose size is well compatible with those of inner bars. This would 
then be the secondary bar of the Galactic double bar system, while the pri- 
mary would be constituted of the COBE/DIRBE bar and the Long bar to- 
gether. 
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building blocks of bars, i.e. the periodic orbits, in 3D. They found 
that the third dimension introduced considerable complexity to the 
orbital structure. Whereas in 2 D it is the orbits of the xi family 
that are the backbone of the bar dContopoulos & PapavannopoulosI 
I198QI : lAthanassoula et al. Il 19831) . in 3D we have a tree of 2D and 
3D families bifurcating from the xi JSkokos et al. l2002bh . Each of 
these families has its own horizontal and vertical extent. Since the 
extent of the box/peanut will in general be determined by a different 
family than that determining the length of the bar, it is natural for 
their lengths to be dif ferent. Furthermore, the ratio of the l engths 
predicted in this way jPatsis. Skokos & Athanassoulall2002h . IS m 
good agreement with that measured for the COBE/DIRBE bar and 
the Long bar (0.8). 

This structure has also been seen in a number of sim- 
ulations, where the bar forms very thin and after a while a 
vertical instability develops and creates the boxy/peanut fea- 
ture (e.g. Binney 1981; Combes et al. 1990; Athanassoula 2005; 
iMartinez-Valpuesta. Shlosman & Helleill2006h . Results from one 
such simulation are shown in Fig.[T](see Appendix for the details of 
the simulation). This is given as an illustration and not as a model 
of our Galaxy. In the upper panel we give the side-on vie\\0 of the 
bar component. For clarity, the remaining disc, as well as the halo 
are not displayed. This shows clearly that the thick part of the bar 
(i.e. the boxy/peanut bulge) is less extended than the thin part and 
that the thin part protrudes on either side of it. It also gives an esti- 
mate of the relative vertical thickness of the inner and outer parts, 
although this could vary from one model to another. Note that the 
thin and thick part are parts of the same bar, and do not consti- 
tute two separate components. Furthermore, both orbital theory and 
simulations show that, in the present scenario of box/peanut forma- 
tion, there must be a thin part of the bar extending further than the 
box/peanut. 

In the lower panel of Fig. [T] the bar is viewed from an angle 
much nearer to end-or[fl. We do not claim that either the model, 
or the angle of the line of sight with respect to the bar major 
axis is necessarily the correct one. This projection, nevertheless, 
illustrates roughly how one could be mistaken into considering the 
poxy/peanut feature and the thin outer part of this bar as two sepa- 
rate component^. 

By visualising a simulation from many viewing angles, it is 
possible to realise the geometry of the object, but this is not possi- 
ble for real galaxies, where only one viewing angle is possible for 
each case. Thus, in near face-on galaxies we can clearly see the bar, 
while in near side-on ones we can see the boxy/peanut bulge. There 
are, however, viewing angles which are near edge-on, but not quite, 
and where both the boxy/peanut shape and the outer thin bar are vis- 
ible. The most interesting inclination range is between 60° and 80° . 
Several examples have already been discussed in the l iterature, such 
asNGC 7582 with an inclination angle of about 65° JOuillen et al.1 
Il997h . NGC 4442 (Betto ni & Galletta 1994 , ) at approxim atelv 72°, 
and M31 at about 77° (I Athanassoula & BeatonI l2006h . In these 
cases it is possible to get information on the ratio of lengths of the 
box/peanut and bar. For M31, where this was specifically measured 



In a side-on view the galaxy is viewed edge-on, with the line of sight 
perpendicular to the bar major axis 

^ In the end-on view the galaxy is viewed edge-on with the line of sight 
along the bar major axis 

^ A short movie, showing this bar from several viewing angles can be found 
in lam. oamp . f r/recherche-14/dynamique-des-galaxies/ 
sc lent if ic- result s /peanut s/milkyway /mo vies- 98 /. The 
username and the password are RefereeArticle and peanut, respectively. 



with the help of cuts J Athanassoula & Beatonll200^ . it was found 
to be ^ 0.7, in good agreement with the value found for our Galaxy. 
iLiitticke. Dettmar & PohlenI (|200Q|) made a detailed morphological 
and photometrical study of a sample of 60 edge-on, also using cuts, 
and found that the ratio of boxy/peanut length to bar length depends 
on the specific morphology of boxy/peanut bulge and, therefore, on 
bar strength. For peanuts, this ratio is 0.53 ± 0.08, for clear boxes 
0.63 ± 0.08 and for box-like shapes 0.71 =b 0.1. According to these 
numbers our Galaxy would be more box-like, but it should be kept 
in mind that these statistics are based on very few objects (21 in 
total, for this measurement) and the scatter quite high. 

The above arguments seem to exclude the possibility that the 
Long bar and the COBE/DIRBE bar are the primary and secondary 
bars of a double bar system, since the bar lengths and length ratios 
disagree strongly both with observations of external galaxies and 
with simulations. Double systems with bars of comparable length 
have never been observed either in any external galaxy or in any 
simulation. It would thus be very hazardous to assume that our 
Galaxy is the only one known to have such a feature. On the other 
hand, the alternative that the COBE/DIRBE and the Long bar are 
parts of a single bar is in good agreement with observations, with 
orbital structure and with simulations. It is thus reasonable to favour 
this second alternative. 

Yet one inconsistency could still remain, concerning the po- 
sition angles of the bar(s). The angle between the major axis of 
the COBE/DIRBE bar and the Galactic Center - Sun line has been 
estimated to be roughly in the range ^ 15° to ^ 30°. First 
observations of the cor responding angle for the Long bar give an 
estim ate of ^ 40° ( Hamniersley et al. I |200Q|: Benjamin et al. I 



20051; iLopez-Corredoira et al. II2007I : ICabrera-Lavers et al. II2008I) . 
while more recent work favours angles around 25° to 35° (Za- 
sowski et al 2011, private communication). The uncertainties are 
such that the observations could be in agreement with the so- 
lution checked here, particularly since structures such as spiral 
arm or ring segments could contribute to the observed Long bar 
signal. In particular, we stress that in many A^-body simulations 
ther e is often, within the inner ring a short, leading segment (see 
e.g. lAthanassoula & MisiriotisI |2002| . and Fig. [21). Furthermore, 
in many A^-body simulations and galaxies seen near face-on the 
boxy/peanut bulge is fatter (i.e. more extended perpendicular to the 
bar major axis) than the thin outer part of the bar. A more accurate 
answer would necessitate a detailed comparison of a large number 
of simulations to the observations. It would be useful to do this 
work using the Marseille library of high-resolution barred galaxy 
simulations run by one of us (EA). 

2.3 Inter-comparing the various models 

For each model we consider three cases. In the first case, the non- 
axisymmetric component describes only the COBE/DIRBE bar or 
boxy/peanut bulge located at an angular separation of 20° with re- 
spect to the Galactic Centre - Sun line. We will refer to this set 
of models as Case 1. From the discussion in the previous section 
we expect this to have a weaker bar than it should, since the outer 
parts, i.e. the contribution of the Long bar, have been neglected. 
Nevertheless, we believe it is useful to consider such cases because 
a large number of such models have been applied to our Galaxy for 
varied purposes. Indeed the five models we are considering here 
have been built by other authors for different purposes. 

In the second case we will include the effect of the Long bar. 
Since our approach necessitates the use of an analytic potential, 
which has not been calculated so far for objects as complex as the 
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bar described in Sect. 12.21 we will model the bar as the superpo- 
sition of two bar models, a vertically very thick one which has the 
properties of the COBE/DIRBE bar and represents the boxy/peanut 
bulge, and a very thin one which represents the Long bar. Trials 
with A^-body simulations show that this is a very reasonable ap- 
proximation, no worst than other approximations standardly used 
in such modelling. We still need to decide at what angle to place the 
bar major axis with respect to the Galactic Centre - Sun line. Since 
the measurements of the distance of the Long bar are quite accurate, 
it would be tempting to place it at 40° . Simulations, however, show 
that in many cases the end of the bar is not symmetric, but extends 
considerably further towards the leadi ng side, in a form reminis- 
cent o f a short arm of ring segment fe.g. lAthanassoula & MisiriotisI 
l2002h . This is also seen in a number of external galaxies. If it is the 
case for the MW as well, then observers would be measuring dis- 
tances both from stars in the thin outer parts of the bar and from that 
leading segment and that would make the bar look as if it were at 
a somewhat larger angle than what it really is. The difference will 
not be big, maybe 5 or 10°. For this reason, we have considered 
many values of the angle between the bar major axis, but will use 
for most displays the value of 30°, while discussing what the effect 
of changing this angle is. Since these two bars are part of the same 
object, i.e. they rotate together, they should have the same pattern 
speed. To model the Long thin bar, we will use the same type of bar 
as for the thick COBE/DIRBE bar, but with different values of the 
free parameters chosen according to the available observations (see 
Appendix). We will refer to this set of models as Case 2. 

As already discussed Sect. 12.21 and also in the previous 
paragraph, we believe that the angular separation between the 
COBE/DIRBE bar and the Long bar is an artifact due to the un- 
certainties of the measurements and the existence of the leading 
extension. Nevertheless, in order to follow these observations we 
will also consider the case where the angular separation between 
bars is of 20° (Case 3). Clearly, since we assume that the pattern 
speeds of the two bars are the same, this situation is dynamically 
unstable due to the forces and torques between the two bars. We 
nevertheless discuss it briefly, to follow the observations. 

In Table [T] we give a brief summary of the main characteris- 
tics of each model. For each model (col. 1), we give the type of 
bar included in the potential (col. 2), the Solar radius (col. 3), the 
value of the pattern speed (col. 4), and for each of the three cases, 
the corotation radius or the distance of the equilibrium point Li to 
the Galactic Centre, r^i, and the three values of bar strength, in 
sequence a, Qh and Qt,Li • The latter are defined as follows: 

• a (cols. 6, 10, and 14) is the ratio of the radial force due to the 
bar's potential to that due to the axisymmetric background, evalu- 
ated at the solar radius and along the COBE/DIRBE bar major axis. 
We thus have: 

dr 

and a = Qr{Ro), where $6 and $o denote the bar and the ax- 
isymmetric potential, respect ively, and Rp is the Solar ra dius. This 
quantity is analogous to qr in lAthanassoula et al] (Il983h . 

• A similar quantity can be obtained if we use the tangential, 
rather than the radial bar force: 

Qtir) = ^ (2) 

^ dr 

where $ is the total potential, $o is its axisymmetric part 



and the maximum in the numerator is calculated over all val- 
ues of the azimuthal angle 0. The maximum of Qt{r) over 
all radii shorter than the bar extent is called Qb (cols. 7, 11, 
and 15 in Table [B and is often used to measure the bar 
strength |Buta. BlocF & Knapen 20031 : iButa. Laurikainen & Salo 
12004^: iLauri kainen. Salo & Buta.2004riButa. Vasvlvev. Salo et al. 
^2005; Durbala, Buta, Sulentic et al. l2009l : lManos & Athanassoulal 
I2OIII) . The radius where this maximum is achieved is defined as 

• Qt,Li (cols. 8, 12, and 16) is the value of the tangential force 
at the Lagrange radius or corotation, Qt,Li = Qt{r = vli). This 
indicator was introduced in Paper III and was shown to correlate 
well with morphological features of the galaxy (such as the axial 
ratio of the rings or the pitch angle of spirals). 

Note that the values of a and Qt,Li are evaluated at the outer 
parts of the bar or the disc, whereas Qb gives the tangential force at 
a radius that lies inside the bar. 

In the left panels of Fig. [3] we show Qr(r) (top panel) 
and Qt{r) (bottom panel) for the five models with only the 
COBE/DIRBE bar, i.e. Case 1. In the inlet in each panel, we make 
a zoom of the region around the Solar radius. The gray strip marks 
the range between the minimum (3.6 kpc for PMM04) and the 
maximum (4.4 kpc for MR09) corotation radius. From the curves 
in the left panels of Fig. [3] and also from the value of Qt,Li, 
we note that model FuxOl has the strongest bar around corota- 
tion, followed, in sequence by MR09, PMM04, DehnenOO, and 
finally GFIO. Around the solar radius, both a and Qt agree in 
the sequence, which from strongest to weakest is FuxOl, MR09, 
DehnenOO, PMM04, GFIO. For models MR09, PMM04 and GFIO, 
the region of corotation is displaced slightly outwards by the intro- 
duction of the Long bar, but it is not affected by any change in the 
angular separation between the two bars, as expected. Indeed this 
displacement is due to the mass of the Long bar and is independent 
of its orientation. In models FuxOl and DehnenOO the bar is only 
modelled as a non axisymmetric forcing, thus introducing no ex- 
tra mass and no displacement of corotation outwards. Comparing 
the curves of the middle and right panels of Fig. [3] at the corota- 
tion and solar regions, we note that the sequence of the models as 
a function of bar strength does not change. Note, however, how the 
introduction of the Long bar increases both Qr and Qt, in all mod- 
els and specially in PMM04, GFIO and MR09, which have a bar 
that comes from a density distribution. 

Note that the value of Qb only reflects the behaviour of the 
tangential force in the inner parts of the bar. If we were interested 
in this region, we would use the value of Qb as the measure of bar 
strength, and the sequence would be as follows: PMM04, followed 
by MR09, FuxOl, GFIO, and DehnenOO. In Case 3, the introduc- 
tion of the Long bar does not change the sequence, although when 
the two bars are aligned. Case 2, the sequence, from the highest 
to the lowest bar strength, is MR09, GFIO, PMM04, FuxOl and 
DehnenOO. 

Since we are not interested in the innermost parts of the 
Galaxy, but rather in a region around corotation, we will use Qt,Li 
as the measure of bar strength. 



3 THE INVARIANT MANIFOLDS 

The models presented in the previous section are composed of an 
axisymmetric component and a non-axisymmetric component, the 
latter described by either a single bar (Case 1) or two-bars (Cases 
2 and 3), rotating clockwise with a constant angular velocity, Q^. 
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Figure 3. Radial profile of Qr (upper panels) and Qt (bottom panels) for the models used here. Left column: Case 1 with only the COBE/DIRBE bar; middle 
column: Case 2 with the COBE/DIRBE and Long bars aUgned; right column: Case 3 with the COBE/DIRBE and Long bars at an angular separation of 20° . 
The inlays show the Solar neighbourhood region in better resolution. In all panels, the gray strip marks the corotation region (see text). 



Case 1 Case 2 Case 3 



Model 


Bar 


Ro 






OL 


Qh 






OL 


Qh 






a 


Qh 


Qt,Li 


FuxOl 


Quadrupole 


8. 


51. 


4.3 


0.022 


0.38 


0.08 


4.7 


0.03 


0.45 


0.1 


4.7 


0.03 


0.43 


0.1 


MR09 


Ferrers 


7.1 


54. 


4.4 


0.02 


0.49 


0.08 


4.7 


0.02 


0.83 


0.09 


4.7 


0.02 


0.50 


0.09 


PMM04 


Composite 


8.5 


60. 


3.6 


0.005 


0.53 


0.08 


3.9 


0.01 


0.61 


0.12 


3.9 


0.008 


0.54 


0.13 


DehnenOO 


Quadrupole 


8. 


51. 


4.3 


0.01 


0.17 


0.04 


4.5 


0.016 


0.23 


0.06 


4.5 


0.015 


0.21 


0.05 


GFIO 


Ferrers 


8.5 


56. 


4.2 


0.003 


0.22 


0.02 


4.5 


0.006 


0.63 


0.03 


4.4 


0.004 


0.36 


0.02 



Table 1. Characteristics of the five analytical models. In the first and second columns, we write the name of the model and the type of bar, respectively. In the 
third column, give the value of the Solar radius and in the fourth column, we give the value of the pattern speed. Then for each case, namely Case 1 (potential 
with only the COBE/DIRBE bar). Case 2 (potential with the COBE/DIRBE and Long bar aUgned), and Case 3 (potential with the COBE/DIRBE and Long 
bar at 20° of angular separation), we give the corotation radius (tl^ ), the radial force at the Solar position (a), the maximum of the tangential force {Qh) and 
the tangential force at the corotation radius (Qt,Li)- The units are kpc for distance and km kpc~^ for the pattern speed. 



In Cases 2 and 3, both bars are assumed to have the same pattern 
speed. We will work in the reference frame where the bar is at rest 
and we use the convention that in this frame, the COBE/DIRBE 
bar is along the x-axis. We concentrate on the motion on the z = 
plane (the Galactic equatorial plane). A ful l 3D study, albeit for a 
simple logarithmic potential, was treated in lRomero-Gomez et al. I 
(l2009h and revealed that the motion in the vertical direction can be 
essentially described by an uncoupled harmonic oscillator, whose 
amplitude is relatively small and that the 3D structures do not affect 
the motion in the z = plane. 

Our theory is largely based on the dynamics of the Lagrangian 
points Li and L2 of a two-dimensional galaxy system. These are 
located where the first derivatives of the effective potential van- 
ish, along the bar semi-ma jor axis, and are unstable saddle points 
teinnev & Tremainell2008h . Each of them is surrounded by a fam- 



ily of periodic orbits, called Lyapunov orbits (lLvapunovl[T94^ . 
Since these orbits are unstable they cannot trap around them quasi- 
periodic orbits of the same energy, so that any orbit in their immedi- 
ate vicinity (in phase space) will have to escape the neighbourhood 
of the corresponding Lagrangian point. Not all departure directions 
are, however, possible. The direction in which the orbit escapes is 
set by what we call the invariant manifolds. These can be thought 
of as tubes tha t guide the niotion of particles of the s ame energy as 
the manifolds (iKoon et al. l200Ql : lG6mez et al. l2004h . In Fig.gl we 
show that from each Lyapunov orbit (light gray thin curve roughly 
in the middle of the panel), emanate four branches: two of them 
inside corotation (inner branches) and two of them outside it (outer 
branches). Along two of these branches (one inner and one outer) 
the mean motion is towards the region of the Lagrangian point (sta- 
ble manifolds), while along the other two it is away from it (un- 
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Figure 4. Dynamics around Li. The four branches of the invariant mani- 
folds (black curves) associated to the Lyapunov orbit (gray curve). The red 
arrows show the sense of circulation, while the black dashed line marks the 
position of the corotation radius. 



stable manifolds). In this theory the stars travel along the orbits 
trapped by the manifolds. 

We need to stress that the terms "stable" and "unstable" do 
not mean that the orbits that follow them are stable and unstable, re- 
spectively. In fact all the orbits that follow the manifold are chaotic, 
but they are in a loose way "confined" by the manifolds, so that they 
stay together, at least for a few pattern rotations, in what could be 
called a bundle. The terms stable and unstable refer to the sense 
of the motion and are related to the saddle behaviour of the equi- 
librium points. We propose that these manifolds and orbits are the 
building blocks of the spirals and rings (Paper II, Paper III). These 
manifolds do not exist for all values of the energy, but only for 
energies for which the corresponding Lyapunov periodic orbit is 
unstable. This means energies within a range starting from the en- 
ergy of the Li or L2 (^j,Li) a nd extending over a region whose 
extent depends on the potential JSkokos et al. Il2002ah . 

The morphological and kinematical characteristics of the man- 
ifolds in external barred galaxies are studied in papers III-V. Note 
that in these papers, the non-axisymmetric component consists of 
a single bar. The main results relevant to this study can be sum- 
marised as follows: 

• The large-scale structure of the galaxy is related to the strength 
and the pattern speed of the bar (Paper III). We find that weak non- 
axisymmetric perturbations produce manifolds of rRi ring shape 
(an inner ring elongated around the bar and an outer ring whose 
principal axis is perpendicular to that of the bar), while strong non- 
axisymmetric perturbations produce spiral arms or other types of 
rings. 

• The range for the shape and size of both the inner and outer 
ring given by the models agree with that from observations, where 
the axial ratio of the inner ring spreads uniformly within the range 
0.6 — 0.95, while the axial ratio of the outer ring falls within the 
range of 0.7 — 1 (Buta 1986). Furthermore, we find a strong anti- 
correlation between the axial ratio of the rings and the bar strength. 

• The default number of spiral arms given by the manifolds is 
two, since they are associated with the number of saddle points 
of the model. The typical orientation is trailing, i.e. following the 



unstable branch of the manifolds, and their shape reproduces the 
characteristic arm winding often observed in external galaxies, that 
is, the spiral arms first unwinds and then returns to the bar region. 

• The formation of the manifolds, and therefore of the rings and 
spirals, depends on the existence of the saddle Lagrangian points. 
These appear when we perturb the axisymmetric potential with 
a non-axisymmetric component. We have applied this theory to 
barred galaxies, but the saddle Lagrangian points can be due to 
other non-axisymmetric perturbations, such as spiral perturbations. 

• In potentials with a strong m=4 component of the forcing, 
manifolds can also account for four-armed spirals. 

• We also studied the behaviour of collisional manifolds using 
a simple model in which the particles within a manifold lose en- 
ergy (Papers III and V). We compared collisional and collisionless 
manifolds and found that they have very similar shapes, at least 
for rings and for the first half turn of spirals. The collisional rings, 
however are thinner, i.e. more concentrated, than their collisionless 
counterparts. 

• Bars are known to evolve secularly, becoming longer and 
stronger while slowing down. As a result, the Lagrangian points 
move outwards and the Qt,Li increase. This has important impli- 
cations for manifolds, whose shape and location will change adi- 
abatically with time. Furthermore it will bring material to the La- 
grangian points, replenishing the mass reservoirs that fill the mani- 
folds. In the case of rings, once material is trapped in the manifolds, 
it can stay there indefinitely, if the model is stationary, while any in- 
crease in the bar strength can trap more material into the manifolds. 



4 THE INNER PART OF THE DISC 

Here we study the morphology and kinematics of the models men- 
tioned in Sect.|2] To compare them morphologically, in Sect. l4.1l we 
show the invariant manifolds, "tubes" that guide the motion of the 
stars, of a given energy level. To study the kinematics, in Sect. 14.21 
we use the orbits trapped by the manifolds, since these are the ones 
that really trace the kinematics. In this case, and in order to make 
the plots clear and clean, we will consider the orbits of minimum 
energy, since for higher energies, the orbits overlap and the area 
they outline grows thicker (Paper I). Finally, in Sect. 14.31 we dis- 
cuss how the manifolds reproduce the observables of the inner part 
of the Galaxy, i.e. the region between 2 — 6 kpc, using the up-to- 
date MW potentials. 

4.1 The morphology 

Figure |5] illustrates how the invariant manifolds change as a func- 
tion of the pattern speed and the bar mass / amplitude for two of 
the bar types considered here, namely the composite bar (top panel) 
and the ad-hoc quadrupole bar (bottom panel). The results for the 
Ferrers bar are similar to the composite bar and can be seen by com- 
paring with Fig. 4 of Paper III. Note that only the COBE/DIRBE 
bar is modelled here (Case 1) and that as we increase the pattern 
speed of a model, the Lagrangian radius moves inwards, and the 
manifolds become attached to the bar. To better understand the 
plots, we give the ratio of the Lagrangian radius over the bar scale- 
length. For the composite bar models, the bar scale-length is fixed 
to a = 3.13 kpc (see Appendix) and the Lagrangian radius is var- 
ied, therefore, the ratio r^/tt increases from left to right from 1.47 
to 2. On the other hand, for the Quadrupole bar models, the ratio 
tl/cl is fixed to 1.25, so the bar length-scale is different in each 
panel. In the latter case, the ratio is inside the range determined 
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Figure 6. Invariant manifolds for the five selected bar potentials and Case 1. From left to right: FuxOl, MR09, DehnenOO, PMM04, and GFIO. The thick solid 
red line marks the position and length of the bar. In the case of the ad-hoc potentials, it marks the bar scale-length. The thin dashed red circle marks the solar 
radius adopted by the respective authors. 



by hy drodynamic simulations r^/a = 1.2 ±0.2 (lAthanassoulal 
Il992h . while in the former, it is outside it. We can not, however, 
extend our search to smaller rijji values because this would entail 
too big pattern speed values ([Gerhard 2010). 

In Fig. [6] we plot the invariant manifolds for all Case 1 mod- 
els that we consider here, that is, up-to-date MW models with only 
the COBE/DIRBE bar and we use the set of parameters that the au- 
thors consider to describe best the MW disc potential. In each panel 
we plot the invariant manifolds of a given energy. The morphology 
does not depend on the energy value chosen, provided this is near 
the Lagrangian point energy (Paper I) 0. The panels are ordered 
from left to right, from the strongest bar model to the weakest bar 
model, namely FuxOl, MR09, PMM04, DehnenOO and GFIO, ac- 
cording to the Qt,Li measure of the force (see Table [T]). All cases 
present both an inner and an outer ring and can be classified as 
rRi ring morphologies. The axial ratio of the rings, however, are 
clearly different from one model to another. In columns 2-5 of Ta- 
ble [2 we give the axial ratio and the major diameter of the inner 
and outer rings. DehnenOO and GFIO have an inner ring axial ra- 
tio that falls wel l within the observational range (0.6 — 0.95) for 
the SB galaxies (lButalll986h . In contrast, the inner rings in FuxOl, 
MR09 and PMM04 are too elongated with axial ratios ranging from 
di/Di — 0.43 — 0.5, respectively. In the particular cases of MR09 
and PMM04, whose bar potential is built from a density distribu- 
tion, the locus of the inner manifolds falls well inside the bar el- 
lipsoid, so that, when self-gravity would be considered, the two 
could merge, as discussed in Paper IV. Therefore, we cannot con- 
sider the inner rings of MR09 and PMM04 as proper inner rings 
surrounding the bar and this justifies their low axial ratio. The bar in 
FuxOl is ad-hoc and we cannot confirm this fact. For DehnenOO and 
GFIO, the axial ratio of the inner ring increases with decreasing bar 
strength, from being oval, di/Di — 0.66 for DehnenOO, to more 
circular, di/Di = 0.77 for GFIO, in good agreement with what 
was found in Paper IV. There are also differences in the size. The 
inner ring in FuxOl, MR09 and PMM04 is very elongated, as men- 
tioned above, not reaching 3 kpc along the minor axis. On the other 



^ The energy chosen for each model is: | (Ej^puxOi -Ej,Li )/Ej^l^ I = 
6.6 X 10-4, \{Ej^Mm9 - Ej^L^)/Ej^L^ = 2.5 x 10-^, 

\{Ej^PMMOA - Ej^L^)/Ej^L^\ = 6.5 X 10-4, |(^J,De^nenOO " 

Ej,L^)/Ej,L, \ = 7.2 X 10-3, and \{Ej^gfio - Ej^l^)/Ej^lA = 
7.2 X 10-3. 



hand, DehnenOO and GFIO have a major diameter of = 4.5 kpc 
and Di = 4.2 kpc, respectively. 

The axial ratio of the outer ring in the five mode ls falls wel l 
within the observational ranse (0.7- 1) for SB galaxies (lButall986h 
and correlates with the bar strength as measured by Qt,Li , as found 
for other models in Paper IV. FuxOl, with the strongest bar, has the 
more eccentric outer ring, do /Do — 0.74, while GFIO, with the 
weakest bar, has a more circular ring with do/ Do = 0.86. The ma- 
jor diameter is around Do = 6 kpc. In FuxOl the outer ring reaches 
6.9 kpc along the major diameter, and in MR09, and DehnenOO it 
is around 6 kpc, while GFIO and PMM04 have the smallest outer 
rings with a major diameter of 5.3 kpc and 5.0 kpc, respectively. 

In Case 2, where the Long bar is aligned to the COBE/DIRBE 
bar, the global morphology is still that of an rRi ring, but there is 
an effect on the shape and kinematics of the rings (see Sect. 14.21 ). 
The axial ratio of the outer ring decreases approximately by 5% 
in the mean, compared to that of Case 1, while the axial ratio of 
the inner ring decreases on average by 20% (see columns 6-9 of 
Table 12]). Changing the angle between the bar major axis and the 
direction of the Galactic Centre - Sun line makes no difference to 
this result. In Case 3, the axial ratio of the outer ring, do/ Do, es- 
sentially does not change compared to that of Case 2. The inner 
ring, on the contrary, gets more elongated with a mean decrease of 
the axial ratio, di/Di, of about 30% compared to the values given 
by Case 1, in agreem ent with the elongated 3 -kpc arm obtained by 
iHabing et al. I (l2006h from the distribution of maser stars in the in- 
ner MW. As discussed in Paper IV this can be due to the fact that 
we include orbits that form the outer parts of the bar. The size of 
the global structure is similar to that of the single bar model (see 
columns 10-13 of Table |2]). This means that there is in general an 
agreement with what is observed in external galaxies. 

4.2 The kinematics 

This section is devoted to the analysis of the kinematics pro- 
vided by the manifolds in the selected models. Fig. |7] shows 
the line-of- sight heliocentric velocity given at each galactic lon- 

S ;itude, hereafter (l,v) diagrani, obt ained from CO observations 
Dame, Hartmann & Thaddeusll200lh . The main features found in 
the inner 2 — 6 kpc of the Galaxy are marked with solid and dashed 
lines, namely the near and far 3-kpc arms and the over-density 
crossing the (l,v) at / = 0° and ^; = km (GMR or Galactic 
Molecular Ring). The rest of the features are related to the central 
part of the Galaxy or to the outer arms, so we will not consider them 
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Model 


di/Di 


Case 1 

Di do/Do 


Do 


di/Di 


Case 2 

Di do/Do 


Do 


di/Di 


Case 3 
Di do/Do 


Do 


FuxOl 


0.43 


4.6 


0.74 


6.9 


0.31 


4.7 


0.68 


7.2 


0.31 


4.6 


0.74 


7.1 


MR09 


0.46 


4.3 


0.75 


6.2 


0.40 


4.7 


0.71 


6.8 


0.32 


4.6 


0.78 


6.9 


PMM04 


0.50 


3.6 


0.78 


5.0 


0.32 


3.9 


0.76 


6.6 


0.41 


3.7 


0.78 


5.5 


DehnenOO 


0.66 


4.5 


0.80 


6.1 


0.54 


4.5 


0.76 


6.3 


0.50 


4.5 


0.82 


6.3 


GFIO 


0.77 


4.2 


0.86 


5.3 


0.72 


4.4 


0.81 


5.7 


0.78 


4.1 


0.82 


5.1 



Table 2. Sizes of the inner and outer rings in the three cases considered, from left to right, only the COBE/DIRBE bar (Case 1), COBE/DIRBE and Long bar 
aligned (Case 2) and COBE/DIRBE and Long bar at cf) = 20° (Case 3). For each of the cases and for each of the models, we give, in the first and second 
columns, the axial ratio and the major diameter of the inner ring, di/Di and Di, respectively, (in kpc), while in the third and fourth columns, we give the 
same values for the outer ring, do /Do and Do, respectively. 



in our comparisons. Note also that these features were observed in 
the CO, while we compare them to collisionless manifolds. This, 
nevertheless, is possible since in Papers III and V we showed that 
the shape of collisional and collisionless manifolds does not differ 
much, at least for rings and for the first half turn of spirals (except 
for the width, which is smaller in the collisional cases). Further- 
more, in external galaxies where it is easy to check, we note that 
the shapes of the gaseous and stellar rings are very similar, with the 
notable difference that the rings in the young stars are thinner than 
the ones in the old stars, in good agreement with our calculations. 
Finally, in external galaxies there are few, if any, purely gaseous 
rings with no stellar component. This would anyway be contrived, 
since it would involve a strong concentration of gas which would 
not form stars. We can thus proceed with our comparisons. 

For each model, we compute the orbits trapped by the mani- 
folds with the energy closest to the energy of Li (minimum energy 
for the manifolds to exist), and we show them in the first column 
of Figures [i [9] and [TOI for the COBE/DIRBE bar only (Case 1), 
the two-bar models aligned (Case 2), and the two-bar models with 
and angular separation of 20° (Case 3), respectively. In the middle 
and right panels of the same Figures, we plot the {l,v) diagram sep- 
arately for the inner and outer branches of the rings, respectively. 
The colours mark the different parts of the ring. We also show in 
a thin dotted black line the axisymmetric component of the termi- 
nal velocity curve, this is, the maximum line-of-sight velocity at 
a given longitude. The black solid and dashed lines are extracted 
from Fig.[7]and trace the position of the near and far 3-kpc arm and 
of the GMR, respectively. 

The kinematic features in the (/,^;)-diagram that should be 
characterised and taken into account for the subsequent compari- 
son with observations fSect. l4.3l ) are: first, where and by how much 
the non-axisymmetric component exceeds the circular terminal ve- 
locities; second, the shape of the lines traced on the diagram com- 
pared to the over-densities on the observed one; third, the non-zero 
value of the line-of-sight velocity at / = 0° for these over-densities 
(expansion or contraction in the radial Sun - Centre direction). 

In several models, we observe that the velocities of the orbits 
in the manifolds exceed the circular terminal velocities at low lon- 
gitudes. For example, the inner rings of FuxOl and MR09 in all 
three cases (one or two bars) exceed the terminal velocities at low 
longitudes (see the second column of Figs.[8j[9]and[Tol). This is due 
to the non-axisymmetric motions, but can in some cases be exces- 
sive (see discussion in Sect. 14.31 ). It is also worth mentioning that 
the introduction of the Long bar changes the {l,v) diagram in al- 



Case 1 Case 2 Case 3 

Model vios N vios F vios N vios F vios N vios F 

FuxOl -15. 13. -30. 30. -27. 28. 

MR09 -46. 44. -108. 109. -59. 58. 

PMM04 -51. 50. -94. 95. -42. 42. 

DehnenOO -7. 6. -15. 15. -15. 13. 

GFIO -6. 6. -13. 12. -9. 9. 



Table 3. Line-of-sight velocities with respect to the LSR at / = 0° . N stands 
for the Near 3-kpc arm, while F stands for the Far 3-kpc arm. Units are in 
km . 



most all cases. In Case 2, where the Long bar is aligned with the 
COBE/DIRBE bar, in all models except for GFIO, the velocities of 
the inner ring exceed the circular terminal velocity curve, due to 
the excess of mass introduced. In FuxOl and MR09, however, the 
maximum velocity is well above the terminal curve. In Case 3, it 
increases the velocities in models PMM04, DehnenOO and GFIO 
to values slightly exceeding the terminal velocities at certain longi- 
tudes. 

Table [3] summarises the line-of-sight velocities along the 
Galactic centre - Sun line, i.e. / = 0°. We see that the range of 
variation for this velocity is 6 — 50 km in absolute value, for 
the models in Case 1 . The introduction of the Long bar in Case 2 
increases the line-of-sight velocity at / = 0° in all models and, in 
Case 3, in all models but PMM04, where it slightly decreases. 

The kinematics along the outer ring also provide informa- 
tion about the deviation from the circular velocity curve. Fig. [TT] 
shows the relative deviation of the tangential velocity in the inertial 
frame from the circular velocity for the orbits in all models with the 
COBE/DIRBE bar. The angle is defined as the azimuthal galacto- 
centric angle with origin on the Galactic Centre - Sun line and mea- 
sured clockwise. The circular velocity at a given radius is given by 
the axisymmetric component in the case of FuxOl and DehnenOO 
models and by the axisymmetric component plus the m = com- 
ponent of the bar in MR09, GFIO and PMM04 models. Note that 
in each case, the deviation from a flat rotation curve is less tha n 
20 km s~^, in agreement with Fig. 6 of iKranz. Slvz & Rixl (1200 ih . 
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Figure 8. (/,i;)-diagrams of the five selected models with only the COBE/DIRBE bar, i.e. Case 1: from top to bottom, FuxOl, MR09, PMM04, DehnenOO and 
GFIO model. Left column: orbits in the (x,y)-plane, the colours showing different parts of the rings; Middle and right column: (l,v) diagram of the inner and 
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Figure 5. 2D parameter study of the Composite and Quadrupole bar models 
for Case 1 . On the x-axis we decrease the pattern speed from left to right. 
On the y-axis we increase the bar strength from bottom to top, from a bar 
mass Mb = 7 X 10^ Mq to = 18 x 10^ M© for the Composite 
bar or from a = 0.005 to o; = 0.15 for the Quadrupole bar. The thin 
black circle in each panel marks the Solar radius. Top panel: the composite 
bar model. The panel with gray background corresponds to model PMM04. 
Bottom panel: the quadrupole bar model. The panels with gray background 
correspond to model FuxOl (top: a = 0.02 and = 51.) and DehnenOO 
(bottom: a = 0.01 and = 51.)- 



4.3 Towards a manifold model of the inner MW 

In the previous sections we analysed in an exhaustive way the re- 
sults given by the manifolds for a wide set of models including 
the possibility of modelling only the COBE/DIRBE bar or both the 
COBE/DIRBE bar and the Long bar. We claim that in general the 
manifolds can interpret the observables both morphologically and 
kinematically. We now compare the results with the MW observa- 
tions and analyse in detail what manifold models are more likely to 
reproduce the 3-kpc arm and the GMR. 

The observed characteristics of the 3-kpc arm are two-fold: 
it is elongated along the bar and it has non-zero velocities at 
/ = 0° (— 53 /cms" ^ for the near arm and 56kms~^ for the 



^ 




50 -50 

l(deg) 

Figure 7. The (l,v) diagram obtained from the observations of 

<Dame, Hartmann & Thaddeu sI l200lh and foame & ThaddeusI bOOSh , as 
plotted by ( iRodrfguez-Fernandez & Combeslboosi) . The line-of-sight ve- 
locity is given with respect to the Local Standard of Rest. The solid lines 
trace the position of some remarkable features such as the locus of the spiral 
arms, the 3-kpc arm and the Connecting Arm. The black dashed lines indi- 
cate the contour of the GMR. The solid cir cles are the terminal velocities 
measurements of iFich, Blitz & Star3 ( Il989h using CO, wh ile the triangles 
are the terminal velocities determined from the HI data of iBurton & Lisztl 
ll993h . The boxes mark the position of the Sagittarius, Scutum, Norma 
and Centaurus tangent points, located, respectively, at I ~ 50°, I ~ 30°, 
I - -30°, and / - -50°. 



far arm JPame & ThaddeusI l2008h . although other studies based 
on ammonia and water in absorption give a value of —43 km 
dWirstrom et a/. 1 1201 d) . This sets a tentative estimate of the error 
bar of about 10 km s~^. Among the set of models considered here, 
we can say that when only the COBE/DIRBE bar is considered 
(Case 1) MR09 and PMM04 reproduce this feature (see Table [3l 
whereas FuxOl, DehnenOO and GFIO have absolute values below 
15 km s~^. 

The introduction of the Long bar, either aligned or misaligned, 
in the models makes the inner ring more elongated and increases 
the velocities along / = 0°, suggesting that the models with the two 
bars have a tendency to better reproduce the observed line-of-sight 
velocity of the 3-kpc arm (see Table [3]). There is a big uncertainty 
in the model parameters, including the angular position of the bar 
in Case 2. Increasing the angle from 20° to 40° makes the line-of- 
sight velocities in the region of the 3-kpc arm to increase in absolute 
value, i.e. increases the rectangular shape of the orbits in the (l,v) 
diagram, while it decreases the size of the necklace- shaped feature 
corresponding to the GMR. In the three selected cases, we can see 
that several models can reach the observed line-of-sight velocities, 
and even exceed them, whereas DehnenOO and GFIO clearly do not 
fit the observations. 

From Figs. [1 [Tol and [9] we observe that FuxOl in all three 
cases and MR09 in cases 2 and 3 have velocities that highly ex- 
ceed the terminal velocity curve and, therefore, they do not satisfy 
one of the main constraints of the (/,v)-diagram. This fact can be 
due to two reasons. First, these models have strong bars, i.e. the 
bar perturbation is quite noticeable in the outer parts. The second 
reason is morphological and related to the previous one. The or- 
bits trapped in the inner ring resemble the cuspy xi orbits of the 
bar, responsible for the maximum in the terminal velocity curves 
teinnevg? al Ill99lh . So the inner branches of the orbits are in the 
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Figure 11. Relative deviation of the tangential velocities in the inertial 
frame with respect to the circular velocity of the outer ring as a function 
of the azimuthal galactocentric angle 0. Top panel: Case 1, models with 
only the COBE/DIRBE bar. Middle panel: Case 2, two-bar models aligned. 
Bottom panel: Case 3, two-bar models at an angular separation of 20° . The 
colours represent the same type of models in all panels. 



immediate vicinity of the orbits of the bar, making an increase of 
the velocity at low longitudes. 

All the above discussion suggests that at least some of the 
model parameters may not be appropriate for the MW, and most 
of them may not be optimum. Indeed, they were selected in the pa- 
pers described in the Appendix so as to give optimum results for 
the COBE/DIRBE bar only, and not for the two bars together. We 
have nevertheless followed them here because they are the standard 
reference models in the field. In a future paper we will search for 
the parameter values which are optimum for the COBE/DIRBE and 
Long bar system. 

The spatial locus of the GM R is n ot well estab- 
lished. IClemens. Sanders & Scovilld (Il988h claim that 
the peak emission of the GM R is at ^ 5.5 kpc, while 
iBinnev & Merrifieldl (Il998h . iDame. Hartmann & ThaddeusI 



around ^Rq. IClemens. Sanders & Scovill a 
the GMR is an almost circular ring with a pitch angle of about 
4°, \yhile other authors (jFux 1999: .Dame. Hartmann & ThaddeusI 



I2OOII : iRodriguez-Fernandez & CombesI 2008h favour the possi- 
bility that the molecular ring is composed of several inner spiral 
arm segments rooted in the central bar. Kinematically, though, 
it is somewhat better constrained. It defines an over-density in 
the (l,v) diagram, crossing the plane from ^ 100 km and 
/ ^ 30° to ^ —100 km and / ^ —40°, passing through the 
origin and having a typical necklace shape. Morphologically, the 
outer branches of the manifolds in most of the cases have a size 
appropriate to the GMR, although FuxOl, MR09 and DehnenOO 
have a major diameter closer to the Sun's position than to halfway 
from the Galactic Center. The shape varies from one model to 
another, the models of DehnenOO, PMM04 and GFIO (with only 
one bar or with two bars) having the shape nearest to circular. 
Kinematically, not all models fit the enhancements related to the 
GMR. The manifold branches plotted in red in Figs. |9] and [TOl 
(Cases 2 and 3) fit well inside the dashed lines in all models, 
while the dark blue branches lie inside the dashed lines only for 
models GFIO and PMM04. FuxOl, MR09 and DehnenOO the blue 
line makes a loop that resembles the over-density f ound between 
/ = 30° - 40° (see Fig. 5 of iRathborne et all (l2009h . where 
they plot the (l,v) diagram given by the ^^CO emission of the 
molecular clouds and clumps of the Boston Galactic Ring Survey). 
Note, in addition, that the angular separation between the bars 
does not affect the global shape of the outer manifolds, either 
morphologically or in the (l,v) diagram. 



5 OBTAINING THE OUTER SPIRAL ARMS 

As mentioned in the Introduction, one of the main features 
of the M ilky Way is the prese nce of two stellar massive spi- 
ral arms (iBeniamin et al. I Eo05h . In this section, we discuss the 
COBE/DIRBE bar models focusing on the parameters necessary 
in order for the manifolds to reproduce a two- armed spiral shape. 

In Sect [3] and Paper II we showed that in a barred galaxy po- 
tential, the two parameters that influence most the shape of the in- 
variant manifolds in the outer parts are the bar mass (or amplitude 
of the m=2 component for quadrupole bars) and its pattern speed. 
In subsequent papers, we analysed the effect of the variation of both 
parameters and showed that the resultant gross global morphology 
is not model dependent, and that, for all models, faster and more 
massive bars produce more open spiral arms, while slower and less 
massive bars produce symmetric outer rings (see Fig. |5] and paper 
III). 

In the previous section, we computed the invariant manifolds 
for five selected MW models with the default parameters chosen 
by the authors. In all cases, the global morphology is that of an 
rRi ring whose extension does not reach the solar radius. As seen 
in Fig. \5\ if we want to reproduce the outer spiral arms of the 
Galaxy using orbits confined by the manifolds, we should signif- 
icantly increase the bar mass or strength, and the pattern speed of 
each model. In Fig. 8 of Paper III, we plot the different models con- 
sidered as a function of the strength parameter Qt,Li- We observed 
that the different morphologies are grouped as a function of Qt,Li, 
so if the bar model has Qt.Li < 0.1 we have Ri outer rings, while 
for 0.1 < Qt,Li < 0.2 the model gives R[ pseudo-rings, and for 
Qt,Li > 0.2 it gives spirals and the other type of rings (R2 and 
R1R2). So we need to use model parameters that give a value of 
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Qt,Li of around 0.2, or larger. Note that this value is significantly 
larger than the values given in Table [T] We will see if this is a plau- 
sible option for the MW. 

In the case of quadrupole bars, the minimum values of a 
and Q.h necessary to obtain tightly wound spiral arms are 0.05 
and 50 /cms" ^ kpc~^, respectively (bottom panel of Fig. [5]). A 
model with two open spiral arms that crosses the solar radius can 
be obtained with a ^ 0.1 and ^ 40 km kpc~^ . Even 
though the pattern speed falls within the range found by simulations 
tenglmaier & Gerhardll 19991) , the bar strength is too high. Remem- 
ber that a is the radial force along the bar semi-major axis and at 
the solar radius. It means that with a quadrupole model and a flat 
rotation curve, we need a bar whose strength is about 10% of that of 
the disc at the solar radius. Besides, the Qb parameter of this model 
is 2.0 8 far outside the observed range jBlock. Buta, Knapen et al. I 
I2OO4I) . If we vary the shape of the rotation curve and we make it 
slightly decreasing in the outer parts (/3 = —0.2), we can obtain 
two open spiral arms with a ^ 0.075 and Qb ^ 57 km kpc~^. 
The radial force of the bar in the solar neighbourhood is still quite 
high but it gives Qb = 0.63, implying a strong ba r, but within the 
observed range (iBlock. Buta, Knapen et al. |[2004h . This might in- 
dicate that if we want to reproduce the spiral arms using manifolds, 
it seems necessary to have a decreasi ng rotation curve in the outer 
parts. A recent study conducted by IXue et al. ] (I2OO8I) analysed 
more than 2400 stars from the SDSS survey and concluded that 
it is indeed possible that the rotation curve of the MW decreases 
in the outer parts, by a factor within the range (3 = (—0.2, —0.1), 
comp atible with the suggested model. Other studies (iBrand & Blitzl 
Il993h suggest that the rotation curve of the Milky Way could be 
slightly increasing in the outer parts, meaning that we should need 
stronger bars in order to obtain open spiral arms with manifolds. 

In the case of Composite bars, the force decreases very 
abruptly in the outer parts of the disc, see Fig. [3] In order to ob- 
tain values of Qt,Li of the order of 0.1 or 0.2 we have to increase 
the pattern speed up to at least Qb 75 km kpc~^ for a bar 
mass of 1.4 X 10^° Mq. If, on the other hand, we decrease the 
bar axial ratio, also related to the bar mas, to b/a = 0.25, we 
can obtain values of Qt,Li of the order of 0.1 with a bar mass of 
~ 2.0 X 10^° Mq and a pattern speed of Qb 70 km kpc~^ . 
Nevertheless, the value of Qt,Li is still too low. The case of Ferrers 
bars is analogue to the Composite bar. 

The main conclusion here is that it seems that in any case the 
COBE/DIRBE bar is not strong enough to make two open spiral 
arms that could reach the outer parts of the disc. Even though the 
introduction of the Long bar makes the inner branches of the man- 
ifolds more elongated (Sect. 14.11 ), it does not change the previous 
conclusions. Models with an extra component to the potential, for 
example a spiral forcing, should be explored in the future as a pos- 
sible option to induce spiral arms through manifolds. 



6 SUMMARY AND CONCLUSIONS 

In this paper, we applied the invariant manifolds theory for the first 
time to the Milky Way. We first presented a thorough discussion 
about whether the Long bar and the COBE/DIRBE bar are, respec- 
tively, the primary and the secondary bar of a double bar system 
in our Galaxy. We dismissed this alternative, because their proper- 
ties are in strong disagreement with the properties of double bars 
in external galaxies and with those in simulations. We then consid- 
ered the alternative that the COBE/DIRBE bar and the Long bar 
are simply parts of the same bar, the former being the boxy/peanut 



bulge and the latter the thin outer parts of the bar. We compared 
the morphology thus obtained with that of external barred galaxies, 
with that obtained from orbital structure and with simulations, and 
found very good agreement in all cases. We also proposed some 
reasonable possibilities to explain the fact that observations indi- 
cate a misalignment between the major axes of the COBE/DIRBE 
bar and of the Long bar. 

We then selected five characteristic models in the literature 
that include the COBE/DIRBE bar, and, to account for recent ob- 
servations and the above mentioned possibility, we have also con- 
sidered the case with an additional Long bar. We analyse the mod- 
els, first in terms of forces, and then we compute the manifolds in 
three cases, namely when only the COBE/DIRBE bar is modelled 
(Case 1), when the two bars are aligned at 30° from the Galactic 
centre - Sun line, as observed in a large amount of external galaxies 
(Case 2), and when the two bars have an angular separation of 20° 
according to the observations (Case 3). Note, however, that in the 
models discussed here, both bars rotate at the same pattern speed, 
which makes Case 3 dynamically unstable. 

Regarding the questions formulated in the Introduction, we 
confirm that the observed features of the inner parts of the Galaxy, 
more specifically the 3-kpc arms and the GMR, can be plausibly in- 
terpreted using manifolds. The morphological and kinematic anal- 
ysis shows that the inner and outer rings described by the manifolds 
can well represent the 3-kpc arms and the GMR of the inner parts 
of the Galaxy, respectively. Not all models, however, are suitable to 
describe both features at the same time and in both ways. In gen- 
eral, the GMR is well reproduced by all models in all three cases, 
but the 3-kpc arms depend more on the characteristics of the poten- 
tial. 

We can tentatively conclude that the two bars of the Galaxy, 
that is, the COBE/DIRBE and the Long bar, are necessary to re- 
produce the observables and, in this case, among the five models, 
the ones with less strong bars, namely PMM04 jPichardo e t al. I 
2Q04), DehnenOO (Dehnen 2000a) and GFIO (Gardner & F Wnnl 
I2OIQI) have a better fit. Another conclusion of this work is that the 
value of the angular separation between the bars does introduce 
changes both in the morphology of the rings and in the {l,v) dia- 
gram. 

We also analysed in detail what bar parameters would be nec- 
essary so that the manifolds have a spiral arms morphology con- 
cluding that, for quadrupole bars, a stronger bar is necessary, while 
for the composite and the Ferrers bars, the force decreases too 
abruptly in the outer parts of the disc, and it is difficult to obtain 
open spiral arms by only increasing the bar mass or its pattern 
speed. In the case, though, where the global morphology would be 
one with two spiral arms, the rRi configuration would be lost. This 
suggests that the most probable solution in the manifold framework 
would be a more complex potential with one or two bars in the inner 
part and a spiral further out (in preparation). 

In a future paper we will revisit the most successful of the 
models considered here and vary their main parameters - i.e. the 
mass and the pattern speed of the bar, and the angle of the bar major 
axis with the Galactic center - Sun line - searching for the values 
that give the best fits to the observations. 
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APPENDIX 

Here we describe the models w e used throughout the paper , namely 
models wi th a quadrupole bar (lFux|[200ll : lDehnenll2000ah. a com- 
posite ba r (Pichardo et al 2004*), and a Ferrers bar ( Ferrers|[l877l : 
iMelnik & Rautiainen 2009; Gar dner & Flvnnll201Ql) . We also give 
a brief description of the simulation used in Sect. 12.21 



The quadrupole bar model: FuxOl and DehnenOO 

The quadrupole model consists of the superposition of an axisym- 
metric component given by a simple power-l aw rotation curve and 
an m = 2 -type potential for the bar, as in iDehnenI feOOOal) and 
If^ ( l200lh . We refer to these models in the text as DehnenOO and 
FuxOl, respectively. 

The potential corresponding to a power-law rotation curve is: 



MR) 



2 
^0 



{2^)-\R/Ror^ 
\n{R/Ro), 



/3 = 0, 



(3) 



where Ro = 8 kpc denotes the galactocentric distance to the Sun 
and vo the local circular speed. Since in this paper we compare the 
two Galaxy models, we will keep the original values of Ro and 
Vo, namely 8 kpc and 220A:ms~^ for DehnenOO, and 8 kpc and 
200 km for FuxOl. The parameter is related to the shape of 
the rotation curve, with /3 = for a flat rotation curve, /3 < for 
a falling rotation curve and /3 > for a rising rotation curve. We 
use flat rotation curves in this paper, unless otherwise stated. In the 
case of falling rotation curves, we use a value of /3 = —0.2. 

The bar is described as the m = 2 component of the Fourier 
decomposition of the potential: ^b{R, 0) — A{R) cos (2^), where 



A{R) = A 




R^Rb 
R ^ Rb, 



(4) 



where Rb and Ab are the size and the amplitude of the bar. In both 
models, Rb = O-Sr^i and r^i =4.35 kpc. This Lagrangian radius 
together with the rotation curve given by Eq. (|3]), corresponds to a 
pattern speed of = 51 km kpc~^. 

The authors measure the strength parameter of the bar by the 
ratio of the forces due to the bar and to the axisymmetric power-law 
component at the Sun Galactocentric radius on the bar's semi-major 
axis. 



Ab fR^y 

v^^ \Ro) 



(5) 



Note that this measure is dimensionless, it is directly related to the 
bar amplitude and it is related to Eq. Q by a = Qr{Ro)- When 
the Long bar is introduced in the potential, it is described also with 
an m = 2 component of the Fourier decomposition of the potential 
with Rbi = 0.92rLi , r^i = 4.35 kpc and ai = 0.6a. 



The composite bar model: PMM04 

The composite bar is an analytical model designed to fit the den- 
sity profile of the bar given by COBE/DIRBE. Again, it consists of 
the superposition of an axisymmetric component and a bar. The ax- 
isymmetric component is the result of the superposition of a bulge, 
a disc and a dark matter halo. The po t ential we used is based on the 
one considered bv I Allen & SantillanI (Il99lh to fit the axisymmetric 
component of the MW. The bulge and the disc are mod elled using 
a Miyamoto-Nagai potential (iMivamoto & Nagailll975h : 



^bi{R,z) 



Ml 



+ ^2 + 52 



,2^1/2' 



(6) 



where (R, z) are the cylindrical coordinates. Mi = 1.4 x 10^° M(r 
is the bulge mass and hi = 0.3873 is the bulge scale-length. 



^d{R, z) = 



M2 



i?2+ a2 + {z^ + hl) 



1/2 



1/2' 



(7) 



where M2 = 8.56 x 10^° M© is the disc mass, a2 = 5.3178 and 
62 = 0.25 are the radial and vertical scale-lengths, respectively. 
The dark matter halo is described using a spherical potential: 



Mr) 



M{r) 
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1.02 



+ In 1 + 



(8) 



where the halo radius is 100 kpc, M3 = 10.7 x 10^° M©, which 
makes the total mass of the halo M(100 kpc) = 8.002 x 10^^ M©, 
and as = 12 is its scale-length. 

The parameters are chosen so that the total mass of the ax- 
isymmetric component is 9 x 10^^ M©, the rotation curve flattens 
at approximately 200 km s~^, setting the galactocentric distance to 
the Sun Ro = 8.5 kpc and the circular velocity at the Sun's position 
to = 220 A;ms"\ 

The bar component is taken from lPichardo et al\ (|2004|) . The 
density distribution is obtained to match the observations from 
COBE/DIRBE. The bar is the result of the superposition of pro- 
late spheroids with density: 



pb{Rs) = po 



sech^{Rs 



sech^ {Rs)e 



Rs ^ Rends 
(I^s-Rer^dc^y \ ,(9) 

Rs ^ Rends 



/ 2 2 2 \ -"-/^ 

where Rs = ( ^ + ^ ) • The parameters ax and ay are the 

ya^ ay J 

scale-lengths of the bar and are fixed to 1.7 kpc and 0.54 kpc, re- 
spectively. The constants Rends hends are defined as abar/dx 
and hend/cLx, respectively. The authors fix abar = 3.13 kpc and 
hend = 0.46 kpc, as the values of the length of the bar and its 
scale-height. The bar is divided in three regions to better describe 
the density. The first two have a fall of sech^ and the third has a 
Gaussian fall that starts where the bar ends. This implies a steep but 
smooth decrease in the density in the outer parts. The bar mass is 
fixed to Mb = 10^° Mq and Qb = 60 km kpc~^ . The model 
parameters of the Long bar, when included in the potential, are 
a = 4.5 kpc, b/a = 0.24 and Mm = 0.6Mb = 6 x 10^ Mq. 
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The Ferrers bar model: MR09 and GFIO 

The third bar model we use is a Ferrers bar as used by 
iMelnik & RautiainenI (l2009h : lGardner & Flvnnl bOlOh . The poten- 
tial is described by the superposition of an axisymmetric plus a 
bar-like component. The axisymmetric component considered by 
the two pairs of authors is slightly different even though both have 

a flat r otation curve. 

In lMelnik & RautiainenI (l2009h , the axisymmetric part of the 
potential has two compone nts: a bulge and a halo. The bulge po- 
tential is a Plummer sphere jBinnev & Trema ine 2008). The bulge 
mass is fixed to 1.22 x 10^° x M© and its scale-length to 0.31 kpc. 
The halo is described by its rotation curve as 



(10) 



where Vmax — 251.6 km is the asymptotic maximum of the 
halo ro tation curve and Vc = S. kpc is its core radius. 

In [Gardner & Flvnnl (l2010h . the axisymmetric part of the po- 
tential has three components: a bulge, a disc and a dark halo. The 
bulge is a superposition of two Plummer spheres, with masses 
3 X 10^ Mq and 1.6 x 10^° M©, respectively, and scale-lengths 
2.7 kpc and 0.42 kpc, respectively. The halo is described by an ax- 
isymmetric logarithmic potential with asymptotic velocity fixed to 
220 km and core radius 8.5 kpc. Finally t he disc results of the 
super position of three Miyamoto-Nagai discs jMivamoto & Nagal 
Il975h with masses: 7.704 x 10^° M©, -6.848 x 10^° M© and 
2.675 X 10^° M©, respectively, and radial scale-lengths: 5.81 kpc, 
17.43 kpc and 34.84 kpc, respectively. The vertical scale-length is 
the same for the three discs and it is fixed to 0.3 kpc. 

As mentioned above, in b oth models, th e bar component is 
described by a Ferrers ellipsoid jFerrersll 18771) , whose density dis- 
tribution is described by the expression 



Po(l 
0, 



m ^ 1 
m > 1, 



(11) 



where m^ — j j}? . The parameter n measures the degree 
of concentration of the bar and po measures its central concentra- 
tion. It is related to the bar mass via the expression: 



r)2n+3 



ixah^p^V(n + l)r(n + 2) /r(2n + 4). 



(12) 



In iMelnik & RautiainenI (l2009l) . the authors fix n = 1, a = 
3.82 A:pc, h = \.2kpc and Mb = 1.82 x 10^° M©. The de- 
fault value for the patte rn speed is Vth — 53 km kpc~^ . In 
iGardner & Flvnnl feOlOh . the authors fix n = 2, a = 3.5 kpc, 
b = 1.4 kpc and Mb = 10^° M©, and the default value for the 
pattern speed is = 55.9 km kpc~^. The model parameters 
for the Long bar, when included in the potential, are a = 4.5 kpc, 
b/a = 0.15 kpc and Mb = 1.1 x 10^° M© for MR09 and 
a 3.9 kpc, b/a 0.15 kpc and Mbi - 0.6Mb = 6 x 10^ M© 
for GFIO. 



The simulation in Sectl2l2l 

The simulation discussed in Sect. I2.2| is part of the library of A^- 
body simulations run by EA, but was not made specifically to 
model the Milky Way, and is only intended for illustration pur- 
poses. In the initial conditions the disc has an exponential horizon- 
tal profile a nd a gec/i^ vertica l one. The halo is described by equa- 
tion (2.2) o f iHernauistI ( [l993h and the bulge by a Hernquist sphere 



(lHernquistlll990l). It was run usi ng the public version of the gyrfal- 
con code (lDehnenll2000blJ2002h . The re solution for Fig jH was en- 
hanced using the technique described by lAthanassoulal (l2005l) and 
made use of the glnemo2 display software, written by J.-C. Lam- 
bert. The same software was also used to create the short movies. 
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